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Abstract 

We present a new model of homogeneous aggregation that contains the essential physical ideas of 
the classical predecessors, the Becker-Doring and Lifshitz-Slyovoz models. These classical models, 
which give different predictions, are asymptotic limits of the new model at small (BD) and large 
(LS) cluster sizes. Since the new theory is valid for large and small clusters, it allows for a 
complete description of the nucleation process; one that can predict the creation of super-critical 
clusters at the Zeldovich nucleation rate, and the diffusion limited growth of large clusters during 
coarsening. By retaining the physically valid ingredients from both models, we explain the seeming 
incompatibilities and arbitrary choices of the classical models. 
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I. INTRODUCTION 



Nucleation refers to the aggregation of identical particles (monomers) into clusters. Its 

n □ 

universality throughout physics, chemistry and biology is well known. References |2J ; 



3|j [llj provide a lineup of the 'usual' (and some unusual) suspects. 



0, 0, 0, 0, 

Also well known are the long-standing challenges that aggre gati on poses to modeling. Two 
classical models of aggregation due to Becker-Doring (BD) ll|, and Lifshitz-Slyozov (LS) 
are incomplete and mutually inconsistent. 

In BD, clusters exchange particles with the surrounding monomer bath by a 'surface 
reaction', and it is assumed that the monomer bath around the clusters has uniform concen- 
tration. This is only possible with infinite diffusivity of monomers. While this description 
is asymptotically accurate for sufficiently small clusters, the uptake of monomers by large 
clusters is strongly controlled by the diffusivity. LS describes cluster growth and shrink- 
age controlled by diffusion of monomers. In LS, the monomer concentration at the surface 
of a cluster is a prescribed function of the local curvature, generally different from the 
'background' concentration, far from clusters. Hence, monomer concentration about a large 
cluster is nonuniform, and there is diffusive transport of monomer into or away from the 
cluster. This physics of LS leads to a prediction for cluster growth that disagrees with BD. 
Furthermore, LS is 'incomplete', in that it does not describe the initial creation of clusters 
from pure monomer. While it is generally accepted that BD is a model for small clusters, 
and LS for large, several questions remain. How to interpolate between the two models? 
What is the characteristic size that separates 'large' from 'small'? What physics governs 
the growth in the intermediate scale? What globally valid model encompasses the whole 
evolution of clusters, from an initial state of pure monomer to the asymptotic self-similar 
distribution of large cluster sizes? 

The current paper presents a new model that retains the essential physical ingredients: 
the clusters gain and lose monomers by a surface reaction that depends on the cluster size 
and the monomer concentration seen on the surface. Monomers outside the cluster undergo 
diffusion with finite diffusivity. These ingredients give rise to a free boundary problem for 
the growth of a cluster that contains a new intrinsic cluster size, fc*, in addition to the well 
known critical size, k c . The critical size, k c , separates shrinking clusters (k < k c ) from 
growing ones (k > k c ). The new cluster size, fc*, indicates the importance of diffusion: the 
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new prediction for cluster growth asymptotes to BD for small clusters with k <C and to 
LS for large clusters with k ^> fc*. In the former case, the diffusion effectively equates the 
surface density of monomer with the far-field density, thus, the surface reaction dictates the 
growth. In the latter case, growth is strongly limited by finite diffusivity. Furthermore, the 
new model of cluster growth interpolates between BD and LS for intermediate cluster sizes 
on the order of fc*. 

The smooth interpolation between BD and LS is crucial for a global model of aggregation 
that describes the whole process, from the initial creation of clusters from pure monomer, 
to the late stage growth-attrition process called coarsening. The essential idea is simple: if 
k c <C k*, as expected in most cases, standard BD describes the nucleation of super-critical 
(k > k c ) clusters and their growth while k c < k <C k*. The super-critical clusters rapidly 
grow to sizes k ^> k*, and their subsequent careers are described by LS. 

Mathematically, we model this physics by a continuum approximation of the discrete 
kinetics. The continuum equations constitute a PDE signaling problem for the distribution 
r(k,t) of large (k ^> k*) clusters in the space of (continuous) cluster size k. At the lowest 
order of approximation, the cluster-size distribution satisfies an advection PDE, in which 
;he growth rate (k vs. k) furnishes the advection velocity. The classical Zeldovich formula 
l3 | . which follows from BD, computes the creation rate of super-critical (k > k c ) clusters. 
Since we assume k c <C fc*, the Zeldovich formula gives rise to an effective source boundary 
condition on k = 0. The initial state of pure monomers is expressed by a zero initial 
condition, r(k, 0) = 0. Information about the amount of small (k < k c ) clusters is not 
expressed directly in r(k,t). Instead, using conservation of particles, we express the amount 
of sub-critical clusters using an integral of r(k,t). 

Our theory does not handle nucleation that happens with k c on the order of fc*. For 
this we suspect that a new theory is needed, one that considers the discrete and fluctuating 
nature of the monomer bath, and does not resort to the diffusion equation, which arises 
from mean-field averaging. 

The paper is structured as follows. In section [H] we present a short summary of the 
classical microscopic aggregation theory (BD). We derive rate constants for attachment and 
dissociation of monomers from a cluster by using free energy and detailed balance arguments. 
The only difference from the classical theory is that it is based on the surface density of 
monomers, the density of monomers just outside the cluster, and does not assume that the 
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monomer density is homogeneous. 

In section 3 we take into account the finite spatial diffusion of monomers. While still 
focusing on a single cluster, we connect the surface monomer density with a far-field monomer 
density, the nearly uniform concentration of monomers far from any cluster. This prescribes 
the growth rate of a cluster as a function of its size and the far-field monomer density. The 
standard assumption in diffusion limited aggregation is that the surface density corresponds 
to a critical cluster i.e., growing and shrinking are equally likely. This seems paradoxical for 
two reasons. First, the free energy of a cluster as a function of cluster size, k, has its global 
maximum at k — k c . On the face of it, this seems to be an unstable equilibrium, but yet it 
is claimed that the cluster remains at the top of this equilibrium. In addition, if that is the 
value of the monomer density, how does the cluster grow or shrink? 

This paradox is another artifact of assuming uniform monomer concentration as in BD. 
It is deconstructed at the end of section 3, by an asymptotic analysis which exposes the 
stabilizing role of finite monomer diffusion: If monomer concentration at the surface of 
a large (k ^> fc*) cluster has large deviation from the critical value described above, the 
surface reactions rapidly absorb or expel monomers. Consequently, due to finite diffusivity, 
the surface monomer concentration undergoes a collateral adjustment towards the critical 
value and the rapid reactions are turned off. The growth rate due to diffusion is much slower 
and dictates the evolution of large clusters. 

In section IIVI we turn to the ensemble of all clusters. If the density of monomers is 
below a certain saturation value, an equilibrium exists, in which large clusters are extremely 
unlikely. For a 'super-saturated' ensemble, with monomer density exceeding the saturation 
value, clusters greater than a 'critical cluster size', k c , have a strong tendency to persist and 
grow. In this super-saturated case, there is no equilibrium; the distribution of clusters is 
continuously changing. Initially, there is nucleation, which is the creation of super-critical 
clusters. The calculation of the nucleation rate based on BD is reviewed here. 

As stated before, this paper proposes a PDE signaling problem for the distribution of 
cluster sizes that quantifies the complete evolution of the aggregation process. Section 5 
contains the assembly of the signaling problem from the component parts in sections 2-4. 
It has a peculiar nonlinearity, in which the advection velocity in the PDE and the boundary 
condition (BC) at k = depend on the monomer density as a parameter. The monomer 
density can be written as an integral of the solution, and herein lies the nonlinearity. 
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The nonlinearity makes the task of solving the equations difficult enough to warrant 
placing it in a separate, forthcoming paper. 



II. CLASSICAL BECKER-DORING MODEL 

Becker-Ddring theory (BD) imposes simplifying assumptions at the outset: The clusters 
are assumed to be uniformly distributed in a dilute 'bath' of monomers. This assumption is 
adjusted regarding the distribution of monomers in the next section, when we add diffusion. 
The clusters are assumed to change size only by losing or gaining one monomer at a time. 
Two large clusters do not fuse together nor does one cluster break into two. This can 
be justified heuristically by noticing that the density of the large clusters is much smaller 
than the (already small) density of monomers, thus the probability of two large clusters 
interacting is small. In addition, the mobility of the large clusters is much smaller than that 
of the monomers, so they are even less likely to stumble upon one another. Similarly, since 
the large clusters have a low mobility (relative to monomers) a cluster that breaks into two 
will, most likely, reconnect quickly, as the two parts remain close together. 

Another important assumption is that the only governing parameter of a cluster is its 
size. The shape of the cluster is assumed to be fixed. This assumption can be weakened to 
require that clusters of same size have the same binding energy and the same surface area. 

To derive the kinetic model of nucleation we introduce the essential quantitative ingredi- 
ents: energy, free energy, and the rate constants of transitions between configurations. 

A. Energy and Free Energy 
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FIG. 1: A schematic cubic crystal surrounded by the monomer bath. 
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The energy of cubic cluster with k monomers is 

e k ~ 3e(k - fc§). (1) 

Here e is the binding energy of a single bond between two adjacent particles, the energy 
needed to break it. 

While the clusters are not assumed to be simple cubes, the general structure of the binding 
energy is expected to remain. Thus, we assume that there is a bulk energy constant, a > 0, 
and a surface energy constant, a > 0, such that the energy of a /c-cluster is 

3 2 

e k ~ k B T(ak - -ah*), for k » 1, (2) 

where, k B T is the Boltzmann factor. The factor of | is added in hind-foresight, as it makes 
some formulas cleaner. Equation @ is only true asymptotically for large clusters. For small 
clusters, the separation between 'bulk' and 'surface' is artificial, and we do not expect (J2J) 
to be quantitatively accurate. In particular, e± = since the binding energy is the change 
in the energy from the unbound state and a cluster with one particle is unbound. 

Next we consider the free energy costs to create a fc-cluster from the monomer bath. The 
bath is characterized by the density p\ of monomer, measured in units of -, where v is the 
volume of each monomer. In other words, p\ is the volume fraction occupied by monomers 

The free energy cost to create a /c-cluster from the monomer bath is 

gk = Sk - k B Tk\ogpi 

~ -k B T Ua + \og Pl )k - ^k'^j , for k > 1. (3) 

Here, k B T 'log p\ is the chemical potential of a monomer in the bath. Rewriting the k ^> 1 
asymptotic form of the free energy gives insight into the existence of a critical monomer 
density, p s = e~ a . We call p s the saturation density of monomers. Setting a = log y in (jSJ) 
gives 

g k ~ k B T ( \ak^> - k log — 
\2 p s 

Thus, when p\ < p s the free energy increases with k, allowing for an equilibrium. When 
Pi > p s , the free energy attains its maximum at the critical cluster size k c , 




log^ 

We investigate the implications of this critical value later. 
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k c ~ I TZZpT I > as Pi -> Pt- ( 4 ) 



B. Kinetics and Detailed Balance 



o 


o 
o 


o 





















J *• 

k 


o 


o 









8 




c 



FIG. 2: The exchange of particles between a cluster and the monomer bath. CkPi is the rate at 
which a monomer gets added to the cluster, and dk is the rate at which monomers leave. 

To study the rate of change in the cluster size we need to model the allowed reactions. BD 
allows two types of reactions: (1) A monomer in the bath can join the cluster; (2) A particle 
on the surface of the cluster can dissociate from it and enter the bath. Mathematically, the 
two reactions are modeled as independent Poisson processes. Let us call the rate at which 
particles leave a {k + l)-cluster dk, and the rate at which particles join a /c-cluster CkPi- We 
work under the assumption that the solution is dilute, that is, p\ <C 1. It is reasonable to 
expect that in the dilute limit, the growth rate is proportional to p\. Thus, we include p\ 
in the growth rate, CkPi, so that both Ck and dk are independent of p\. 

Initially we assume a uniform pi, so it is a global parameter. When we add diffusion, we 
are a little more careful and take p\ to be the density at the surface of the cluster. 

Since the dissociation happens on the boundary of the cluster, we expect dk to be pro- 
portional to the surface area, which, in turn, is proportional to /c 2 / 3 , 

2 

dk = ojk'i . (5) 

Here, u includes a per-particle dissociation rate and a geometric factor. 

A standard detailed balance argument relates to dk- Let pi be the value of monomer 
density so that a /c-cluster is in equilibrium with the monomer bath. That is, the cluster 
has no net tendency to grow nor shrink, hence, the adsorption rate, piCk, should match the 
emission rate, dk- The physical requirement for the equilibrium is that the free energy is 
unchanged when a monomer is taken from the bath and added to the cluster, 

9k+i - 9k = + £fc - k B T log pi = 0, 
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hence, 



pi — e k B T 



The detailed balance relation between c k and d k is therefore, 

c k = e k B T d k . (6) 

We use the adsorption and emission rates, c k p\ and d k , to derive a kinetic equation for the 
expected change in size of a cluster. In a small time span 5t, the expected change in cluster 
size, 6k, is: 

(6k) = (c kPl -d k )6t. (7) 

Using the model for the binding energy of large clusters (J2J) and relation (JHJ) between c& and 
<4, equation (J2J) can be rewritten as 

~ |^ Pl ~ Pg A;§ - ak^j 8t, for jfe > 1. 

In the appendix we show that if the variance in k is much smaller than (k) initially, it will 
remain so, as long as (k) is bounded away from the critical size k c in (jlj). In this case we 
approximate the evolution of k(t) for a given cluster as deterministic, and governed by the 
ODE 

k = uj (r)k% - ak^) . (8) 



Here, rj is the super- saturation, defined by 



Pi - Ps , a , 

V = • 9 

Ps 

The super-saturation in Q is generally a function of time, due to the exchange of particles 
between clusters and the monomer bath. By conservation of the total particle density, the 
average value of p\, and hence rj, is determined from the densities of all clusters with k > 2. 
Within the framework of BD, which assumes that the monomer density is uniform, we simply 
set rj to this average value, and, in this sense, (JHJ) is the BD prediction for cluster growth. 
However, if the diffusivity of monomers is finite, the density of monomers seen at the surface 
of the cluster will be different from the average value far away. We propose that (JHJ) holds 
generally, with rj equal to the value of super-saturation seen at the surface of the cluster. 
Equations (JHJ) and expose the kinetic significance of the saturation density p s = e~ a , and 
the critical cluster size k c in (J3J), which we rewrite using rj, 

k c , as rj - + . (10) 



If the surface value of p\ is less than p s , all the clusters shrink, regardless of k. For p\ > p s , 
i.e. rj > 0, the critical size, k c , separates growing, super- critical clusters (k > k c ) from 
shrinking, sub-critical ones (k < k c ). 

It should be noted that while the expected change in size of sub-critical clusters is nega- 
tive, there is a small probability for a sub-critical cluster to grow and become super-critical. 
The rate at which this happens is estimated by the Zeldovich formula, to which we come 
back in section IIV1 First we add finite diffusion of monomers to the model and see how it 
affects the growth rate. 




FIG. 3: The cluster is surrounded by a inhomogeneous monomer bath. The flux of monomers in 
the bath is diffusive, and the cluster reacts to the local monomer density. 



III. ADDING MONOMER DIFFUSION 

Diffusion is the usual model of transport for monomers in the bath. For a finite diffusion 
coefficient, D, the exchange of particles between a cluster and the bath directly outside of it 
leads to non- uniform monomer density. Hence, the monomer concentration, p 1; is a function 
of position and time that satisfies the diffusion PDE. 

To add diffusion to the BD model, we make a couple of additional assumptions. The 
cluster is taken to be spherical, filled with monomers, each taking volume v. Thus, the 
number of particles in the cluster, k, and the radius of the cluster, a, are related by 

kv = fa*. (11) 

The monomer density is assumed to be radially symmetric (with the cluster centered at the 
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origin.) Thus, the density, p\(r,t), satisfies the radially symmetric diffusion PDE in R 3 

d t pi = ^<9 r (r 2 <9 r pi) in r > a, (12) 

There are two BC at r = a: one is the kinetics due to BD as in equation (jHJ). That is, 
pi (a, t) is related to k and k by 



k = w\7i(a,t)h* -ah*), (13) 

where, 

p 1 (a,t)-p s 
r}{a,t) = - 



Ps 

is the super-saturation seen at the surface of the cluster. 

The second BC results from conservation of particles. The particles that are added to 
the cluster can come from two possible sources: particles in the solution surrounding the 
cluster, which join the cluster as it engulfs them, and particles added by the diffusive flux. 
This can be put in a simple equation: 

d(l -pi(a,t)) = D(d r pi)(a,t) = Dp a (d r rj)(a,t). (14) 

Finally, there is a BC at oo: the monomer concentration has the asymptotically uniform 
value poo far from the clusters, 

Pi — > Poo as r —> oo. (15) 

In the context of the full aggregation problem, p^ is generally a function of time, which 
follows from overall conservation of particles. Equations ()llH15j1 constitute a free boundary 
problem (FBP) for a(t) and pi(r, t) in r > a. 

The analysis of the FBP begins by identifying suitable non-dimensional variables. We 
assume that the local super-saturation 

V{r,t) = . 

Ps 

is uniformly small in r > a. Hence, we introduce e as a gauge parameter for r)(r,t) and we 
replace rj(r,t) in (|llfJT5j) by er](r,t). The analysis of the FBP is based on an e — > limit 
process. However, we do not yet take the e — » limit. It remains to determine the scaling 
of the other variables k, r, a and t with e. These follow from simple physical balances: in 
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the 'kinetic' BC ()13|) . r](r, t) has order of magnitude e, so the two terms of the RHS balance 
when k has magnitude (|) 3 . Notice that this is the critical cluster size (fTUj) for r/ = e. The 
scalings of r/ and k are recorded in the scaling table 

Scaling Table 
Variable rj k a,r t 



Unit e ( Z Y -v 1 / 3 

The unit of r and a is the radius of a cluster with k = (7) 3 particles. The unit of t is chosen 
to balance the LHS of the BC (|13|) with the two terms on the RHS. 

The equations of the non-dimensional FBP that follow from (J11H15)) are 

k = ^a 3 , (16) 

9tV = — Tn — \d r (r 2 d r r]), in r > a, (17) 

k = 7 1 (a J t)k 2/3 - k 1/3 , (18) 

D £p s 
ojv 2 / 3 a 

V — ► ^00 as r — > 00. (20) 



afl-p s (l + £77(a,t))) = _-!^(9 r?7 )(a,t) (19) 



Here, 77^ is the asymptotically uniform value of r](r,t) as r — ► 00. The characteristic time 
of its variation is assumed to be comparable to or larger than the unit of time from the 
scaling table. 

The dimensionless constant J^ ri can be interpreted as a ratio of characteristic times for 
two different physical processes. Recall that - is the characteristic time for a monomer on 

2/3 

the surface of a cluster to dissociate into the solution. The ratio is the characteristic 
time for a monomer to diffuse a distance comparable to its own size. The conventional 
assumption is that the "dissociation time" is much longer than the "diffusion time," so that 

1 v 2 / 3 D 

In this limit, the diffusion equation (|17|) reduces to a radial Laplace equation 

d r (r%r)) = 0. 

The solution with 77 = 77^ as r — > 00 is 

ri{r,t) = rj 00 + {rj(a,t)-rj 00 )^-. (21) 
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The time dependence of r\ is implicit due to the time-dependence of its values at r = a and 
r = oo. 

In the dilute limit, with p s <C 1, the conservation equation ()19|) reduces to 

D ep s 



Substituting 77 from (J2TJ), this is becomes 

D sp 

ujv 2 / 3 a 



■(d r rj)(a,t). 



ah = ~ V{a,t) ) • (22) 



Equation can be converted into an equation for k using (|TS|). This gives 

fc = efifyoo- r)(a, t))k*, (23) 
where /x is the dimensionless constant: 

Since /x is a product of a large number, , and a small number, ^, it can be large or 
small. We therefore entertain any value of p. 

The two equations for k, (fTHJ) and (J23J), involve the super-saturation at the surface of the 
cluster, r)(a,t). Solving for k and fj{a,t) gives 

n h 2 / 3 - k 1 / 3 
k = J , (24) 

Voo + — 

»KM) = — ( 25 ) 

1 + — 



ODE (j2H) indicates a second characteristic cluster size besides k c . In units of (|) this 
cluster size is 

fc* = (e/x) 3 . 

In the original variables, A;* is a combination of basic physical constants: 



3 J3 



^ = (a/i) 3 = (3-167r 2 )^-4. (26) 

Notice that for k <^k* equation (J23jl asymptotes to BD. For k ^> fc*, the asymptotic form 
of J2H) is 

fc ~ ejz (rfco^ - l) . (27) 
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Restoring original units, (|27j) becomes 

k = d (rjook* - a) , rf= (3- 16tt 2 )5^, (28) 

VT> a 

which is the standard result for diffusion limited growth (DLG) |3j- Equation (|25[) shows 
how the surface value of super-saturation differs from the uniform value, 7700, far from the 
cluster. Notice that it is a function of k. We convert it into an equation for Sf., the value of 
monomer density seen at the surface of a /c-cluster (again, in original units): 

This will be important, when we examine the whole ensemble of clusters, and formulate 
evolution equations for cluster densities. 



A. Physical Meaning of k* 

We show that k* is the characteristic size of clusters for which finite diffusion induces a 
significant relative difference between rj^ and rj(a, t). That is, 5rj = ri^—ri^a, t) is comparable 
in magnitude to r]^. A simple examination of two physical balances is sufficient. First, the 
cluster's growth rate balances the diffusive influx of monomers. This is expressed by 

k = a 2 D^^. (30) 
a 

Here, the equality means 'order of magnitude balance'. In the RHS, ^-Srj is the difference 
between monomer densities at 00, and on the surface, expressed in the conventional unit 
of 1 /volume. For quasi-static diffusion, the diffusion zone about the cluster of radius a has 
thickness a, so — estimates (d r pi)(a,t) and the influx of monomers per unit area into 
cluster is estimated by multiplying this by D. Finally, multiplying by the area, proportional 
to a 2 , gives the cluster growth rate, k. 

Second, the magnitude of k as dictated by the surface reactions (fT3"| is 

k = u7]k%. (31) 

Enforcing the equivalence of (J30|) and (jHTJ) and using k = va 3 (order of magnitude equality 
again), we find 

Sr] _ uv 2/3 1 _ / k \ 5 
77 Dp s \KJ 
We see that Sr] is comparable to r] when k is comparable to fc*. 
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B. Critique of DLG and its 'paradox' 

We briefly examine the 'traditional' derivation of ODE ()27|) for DLG, within the frame- 
work of the non-dimensional free boundary problem ()16H2(Jj) . Given r](a,t), equation (J23j) 
gives the growth rate of the cluster that follows from diffusive flux of monomers. In the 
traditional analysis of DLG, r](a,t) is chosen so that the cluster is in equilibrium with the 
monomer bath that surrounds it. Under the current non-dimensionalization, this 'critical 
nucleus' BC reads 

r](a,t) = k~k (32) 

Substituting §2$ for rj(a,t) in (J2EJ) leads directly to the ODE (|27jl . 

By inspection, we see that the traditional equations (|2"T| l3~2*|) arise by taking the efi — > 
limit of f)24l |2~5*)) with k fixed. The alternative limit process, — > oo with A;* = ep fixed is 
more physical: the value of e\i is set by material properties and initial conditions, and we 
expect that clusters eventually grow to sizes k 3> A;*. We have already seen that the ODE 
(124)1 for k converges to the DLG result in this limit, but the expression (125)1 for the surface 
value of monomer density does not converge to the DLG boundary condition (|32~|). Instead, 

r)(a,t) ~ (1 + /^Ari, (33) 

which has an additive term fir/oo in the prefactor of A; -1 / 3 not present in (|32j) . A mathematical 
critique of the 'critical nucleus' boundary condition ()32j) is simple; it results from formally 
neglecting k in the LHS of the 'surface kinetics' boundary condition (|18j). In our result, k 
balances the RHS, even in the limit j- — > 00, resulting in (|33j). 

Recall that the 'critical nucleus' boundary condition in traditional DLG looks paradoxical 
because the 'cluster sits on top of a free energy maximum'. A 'lazy' deconstruction might 
say: "Nothing to explain, the critical nucleus boundary condition is simply incorrect in the 

(more physical) limit ^ > 00 with A;* fixed." Another easy explanation looks at the free 

energy. The free energy Q refers to a simple cluster surrounded monomers of uniform 
density, whereas the actual kinetics we consider involves a non-uniform density pi(r,t) in 
r > a. The actual free energy takes into account the functional dependence of pi(r,t) in 
r > a. These remarks indicate that the 'paradox' in its original form is naive. Nevertheless, 
it points to some physics that is not expressed in the quasi-static model (|24l I25j) of cluster 
growth as it stands. 
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Suppose that we place a cluster of size k into a uniform monomer bath that has the 
'wrong' monomer density, not equal to the surface value Sk given in (|29j). In order for our 
model to be plausible, the surface value of p\ should rapidly relax to Sk in (|29|). We now 
show that the full free boundary problem (|lbH2L)j) implies such a relaxation transient. 



C. The Stability of the Free Boundary Problem 

The relaxation transient is characterized by a balance of the time and space derivatives 
in the diffusion PDE ()12j) . Hence the characteristic time of the relaxation transient is 

tr = (34) 

where a is the cluster radius. The relative change of the cluster radius in this characteristic 
time is small: from ()14j). d has the order of magnitude 2£sl_ Hence, the relative change of 
radius in time t r has magnitude ep s . The small relative change in cluster radius means that 
the cluster radius is asymptotically constant during the relaxation transient, and it remains 
to derive from the full free boundary problem ()llflTK|) a reduced boundary value problem for 
77 (r, t) in r > a, with a fixed. We use the previous units in the scaling table for all variables 
except time t. For t we use t r in ()34|) with a replaced by the characteristic cluster radius, 
jV 1 / 3 . The reduced boundary value problem is 

d t r] = —d r (r 2 d r r]) in r > a, (35) 



r 



\{d r r]) (a, t) = r](a, t) - k~ 3 , (36) 
rj(r, t) —> 7700 as r —> 00, (37) 

in the limit e — > 0, and A = (jA 3 e[l fixed. The far-field super- saturation, 77^, is assumed 
to vary on a characteristic time much longer than t r , so 77^ is effectively constant. 

The time-independent solution of (|3*oH3Tj) for r)(r, t) is (j21j) with rj on r = a given by ()25|). 
We show that this time-independent solution is asymptotically stable. We notice that 

E^-J r 2 (d r r ] ) 2 dr + Y ( V (a,t)-k-^ (38) 

is a Lyapunov functional for equations (j33H3*7j) . The time derivative of E, 

POO 

E = — A / r 2 rj 2 dr, 
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is found by time-differentiation of ()38|). integration by parts, and use of the PDE (|35j) and 
BC f!36|) . Since E is positive definite, and E < 0, it follows that r/(r, t) converges to the time- 
independent solution. We conclude that if the surface monomer concentration is initially 
different from the quasi-static value (J25)) . it relaxes to it in characteristic time t r . 



IV. EVOLVING DISTRIBUTION OF CLUSTER SIZES 

The kinetics equation (J24)) requires 77^ , the super-saturation far from any cluster. To find 
rjoo we look at the joint evolution of all the clusters, each assumed to follow the dynamics 
in (j23)- The clusters are coupled by the combined effect they have on the monomer density, 
and consequently, on the super-saturation. 

Let pk{t) be the average spatial density (in units of -) of A;— clusters at time t. Assuming 
that the total particle density has a fixed value, p, the pk satisfy a particle conservation 

equation. Since a fc-cluster is made of k particles, the total particle density, p, must satisfy 

00 

P = J>p fc . (39) 
fc=i 

With the help of (J39|) . the space averaged super-saturation can be rewritten as a function 
of the cluster densities pk with k > 2: 



Pi - Ps P~ Ps 



00 

(40) 



Ps Ps Ps k 

In the dilute limit with inter-cluster distances much greater than cluster radii, we expect 
that the super-saturation is asymptotically uniform, throughout most of the monomer bath 
far from clusters. In this case, that asymptotically uniform value, ?7oo(^), should be well 
approximated by the spatial average (j4*Uj) . 

Voo = - — — y~] kp k . (41) 

We turn to the evolution of the densities. The pk obey kinetic equations associated with 
the reactions 

/c-cluster + monomer (k + l)-cluster. 

The equations are 

Pk = jk-i ~ jk, (42) 
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for k > 2, where the discrete flux jk is the net rate of creation of a (k + l)-cluster from a 
fc-cluster, 

jfc = c k s k p k - d k p k+ i. (43) 

As before, <ifc is the rate constant for shedding a monomer from the surface of a (k + 1)- 
cluster. For fc 3> 1 it has the asymptotic behavior (jSJ), proportional to surface area. The 
prefactor c k Sk of p^ in (J43|) is the rate constant for adding a monomer. Recall that s& is the 
value of monomer density seen at the surface of a fc-cluster, given by (|2*H|) . and c?; is related 
to e4 by the detailed balance condition (JHJ). A more explicit formula for jk displaying the 

/c-dependence of surface monomer concentration and detailed balance is 

( £ fc+i- £ fc \ 
e k B T s k p k - Pk+i j • (44) 

Equations (|421 14^j) can be summarized as discrete advection-diffusion equations, 



Pk + D- 



£ k+l- £ k 

d k [ 1 - e k B T s k ) p k - d k D + p k 



for k > 2. (45) 



Here D + ,D~ are, respectively, the forward and backward difference operators. In these 
equations, the surface monomer density s& contains the super-saturation 77^ as a parameter, 
and 7/00 is connected to the pk according to f)41jl . So we see explicitly how the cluster densities 
are coupled to each other via the super-saturation. The /c-dependence of s k induced by the 
finite diffusivity of monomers is the essential difference from classical BD. We recover 
classical BD by taking p — > 00, which in turn results from ^ 3 — > 00. Then Sk in (|2T?|) 
reduces to p s (l + 7700), which is the uniform value of pi assumed in classical BD. 



A. Equilibrium 

Equilibria are time independent densities, pk, so that all the fluxes jk are zero, and the 
sum (I39J) . which gives the total particle density, is convergent. Setting jk = in (j44j) gives 
a recursion relation that determines pk from p 1; 

p k = p\e^ , for fc > 2. (46) 

Here, we used Sk = Pi for all k, since pi should be uniform in the equilibrium case. Substi- 
tuting these pk into (j2HJ) gives, 

00 

p = p(pi) = E^i e *- ( 47 ) 
fc=i 
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Hence, equilibria exist for monomer densities p\ so that this series converges. For large 
values of k, the binding energy can be written as 

e k ~ k B T (ah - l^k 2 ^ 

Therefore, convergence happens for 

pi < e~ a = p s . 

In other words, equilibria exist only if the super-saturation is non-positive. The largest value 
of total particle density for which there is equilibrium is obtained by setting p\ = p s in (|47|). 
We denote this critical particle density by p c , 

Pc = p{p s ). (48) 
Since p s <C 1, the first few terms of the series give a close approximation of p c . 



B. Zeldovich Nucleation Rate 



For positive super-saturation 77, there is the critical cluster size k = k c , for which the free 
energy cost to assemble a /c-cluster from dissociated monomers is maximized. In the small 
super-saturation limit 77 — > + , and k 3> 1, it follows from (J3J) and Q, that 

gk~k B T[ %aki -acq] , (49) 



2 

and that the free energy cost (in units of k B T) of the critical cluster is asymptotic to 

For small super-saturation, the free energy cost is high, and an analogy with the famous 
Arrhenius rate suggests that super-critical nuclei with k > k c are produced at a rate pro- 
portional to the exponential e~ 9 . Since this proposed creation rate is exponentially small as 
77 — > 0, one might expect that after some initial transient, quasi-static but non-equilibrium 
values of p k are established for k on the order of k c , in which the discrete fluxes jk in (J43|) are 
asymptotically equal to a uniform value j. This j is the creation rate of super-critical nuclei, 
proportional to e~ 9 . These essential ideas of nucleation kinetics are set forth in a famous 



paper of Zeldovich, on the nucleation of vapor bubbles for under pressurized liquid 
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121. His 



starting point is a discrete system of kinetic ODE's like BD, but he first passes to a PDE 
limit of the ODE's and calculates the nucleation rate from the PDE. Here, we implement 
the essential Zeldovich ideas, but within the framework of the discrete BD ODE's. 

We work in the limit k c <C fc*, so for k on the order of k c there is negligible difference 
between the surface value, s^, of monomer density in (|29|). and the uniform value, pi, far 
from clusters. We show this: in equation (j2HJ) for Sk, we see that Sk ~ p s (l + r/oo) = Pi if 
i < r/oo and ^ < 1. For k c = {o/^f in (HDD and K = (<xp) 3 in (J2SD, we find £ = 
So, tt <C 1 implies - <C r/oo. Furthermore, for k on the order of k c , ^-^ is on the order of 

1/3 

< 1. 



fee 

fe» 

In (J44J) . we replace by pi, and by j, to obtain a recursion equation that determines 
the pk for k > 2 from pi. We write it as 

p±-p^=L ( m. ( 5o) 

Here, gu is the free energy cost of a fc-cluster, given in (J4*9~j) . and p^ denotes the solution pUj) 
of the homogeneous recursion relation with j = 0, and p\ = p\. For positive super-saturation 
?7, pfc — > oo as A; — > oo, and we expect that the p/c in (f5U|) have |£ — > as k — > oo. Summing 
(|5()|l over gives a formula for j. On the LHS, we get a telescoping sum with value 

lim^ = l-0 = l. 

Pi Pfc 

Hence, 

l=jS^—e&. (51) 

In the RHS, g% decreases linearly with k as k — > oo, so the series converges. In the small 
super-saturation limit 77 — > + , we expect that the sum on the RHS is dominated by terms 
with k near k c ~ , where g^ attains its maximum. The relevant approximation to g^ 
as 77 — > + and /c is on the order of fc c is given by (|49|). Also, cZ^. ~ cuk 2 ^ 3 as in (J3J). Hence, 
((nij) has the asymptotic approximation 

l = lj2 k ~~ 3 ^ ak *~ kV - (52) 

^ fc=2 

The final step is the approximation of the sum by an integral, and evaluation of the 77 — > + 
limit by the saddle point method. This leads to the approximation of j, 

3 ~ U \ T- 6 '^ 7 ■ ( 53 ) 
V 07T 
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V. ADVECTION SIGNALING PROBLEM 



We examine the aggregation process, starting from a super-critical density of particles, 
p > p c , all in the form of monomers at time t = 0. That is, 



There is an initial transient, called ignition, in which the first supercritical clusters appear. 
A detailed analysis of the ignition transient appears in a paper by Bonilla et al. Q|. Here 
we give a brief summary. First, sub-critical (k < k c ) clusters are created with quasi-static 
densities close to the values pk in (jJBJ). Of course, the value of p\ becomes less than pi(0) = p, 
since these sub-critical clusters are created from monomers. Hence, the appearance of the 
sub-critical quasi-static densities is accomplished by the decrease of super-saturation from 
an initial value of 77(0) = ^y^- to a smaller value, which we denote by 77*. In appendix IA 21 
we show that 77* is related to p — p c by 



as p — > p+. Here, p = p(pi) is the equilibrium relation (JTTj) between p x and p for < p\ < p s . 
Figure 0] is a visualization of the relation (|54j) between 77* and p — p c . The 'completion' of 
the quasi-equilibrium in 1 < k < k c = (jf^J is accompanied by the appearance of the first 
super-critical clusters with k > k c . The rate of creation rises from zero to the Zeldovich rate 
111 (P| with rj = 77*. We assume 77* is so that k c <C A;*, so the Zeldovich rate indeed applies. 

Now our focus shifts to the evolving distribution of the super- critical clusters. The model 
we present here has three physical ingredients. Two of them are the growth of the clusters, 
and their creation. Both processes contain the super-saturation as a parameter, and the 
remaining ingredient is the connection of super-saturation to the distribution of cluster 
sizes. 

First, growth. We expect a predominance of super-critical clusters with k 3> k* that 
undergo diffusion limited growth. Here is the heuristic argument for not resolving size scales 
comparable to k* or smaller: once a cluster achieves super-critical size, with ^> if 
(see the appendix IA 1|) it continues to grow nearly deterministically, as shown in section 3. 
Since the Zeldovich rate is exponentially small in 77, an exponentially long time elapses 
before the super-saturation shows an significant decrease below the effective initial value, 



Pl {0) = p>p c , p fe (0) = for k > 2. 




(54) 
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Pi(0) 



Pc 




FIG. 4: The darkened curve is the graph of p(pi) in < pi < p s . The dashed line is its linear 
interpolation into p a > p s . It intersects p = pi(0) > p c at pi = p s (l + 77*). The arrow labeled i 
represents the decrease of p\ during ignition, and the arrow a represents the decrease during the 
subsequent aggregation process. 

77*, established during ignition. In this exponentially long time, we expect the super-critical 
clusters to grow to sizes k ^> k*, the regime of diffusion limited growth. 

We assume that the characteristic cluster size k corresponding to significant variations of 
the densities pk for k 3> fc* ^> 1 is itself much larger than fc*, and this motivates a continuum 
limit, 



Here, r(k,t) is a smooth function of its arguments, with the characteristic scale k oik much 
larger than fc*. Substituting for p k into the discrete advection-diffusion equations (Jj3j) 
and using the assumed largeness of k, it follows that r(k,t) asymptotically satisfies the 
advection PDE, 



Here, the advection velocity u = u(k,rj) is identified from ODE for diffusion limited 
growth. We have, 



where rj = rj(t) is the 'background' super-saturation, far from any cluster. 

Next, creation. In the analysis according to Zeldovich, recall that the discrete flux in 



p k (t)~r(k,t). 



(55) 



d t r + d k {ur) = 0. 



(56) 




(57) 
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is asymptotically uniform for k on the order of k c , with value j given by (|o^|). Here, we 
make the stronger assumption that the range of k with the asymptotically uniform value j of 
jk extends to a scale of k much larger than k* but smaller than the characteristic cluster size 
k associated with the continuum limit (J55|) . In this case, we expect an asymptotic matching 
between the continuum limit of jk, given by 

j k ~ u(k,rj)r(k,t), (58) 

and the uniform value, j, of jk, in some overlap domain of cluster sizes k much larger than 
fc*, but much smaller than k. Since k ^> &;* ^> k c = (jj-^j i n the overlap domain, the 
dominant component of u in (j5Tj) is d i] k 1 /^ and (}58|) reduces to 

j k ~ d7]k 1/3 r(k,t). 

in the overlap domain. Hence, we propose the effective boundary condition (BC) 

drjk 1/3 r(k 1 t) — > j = uo^^-e'^n 1 as k — > 0. (59) 

If 7] = rj(t) is known, the advection PDE (JSTIj) and creation BC ()59|) lead to a simple 
determination of r(k,t), starting from the initial condition of pure monomer, r(k, 0) = in 
k > 0. 

It remains to connect r](t) to r(k,t). In the particle conservation identity we now 
distinguish between sub-critical and super-critical clusters sizes. That is, 

p = ^kpk + ^2 k P k - 

The sub-critical sum is approximated by substituting quasi-equilibrium values for pk 
based on super-saturation rj, 

Pk~h = (l + v) k p k s e^, 

and then taking the limit rj — ► 0. The details here are a re-run of the calculation in appendix 
I7T21 We get 

k c 

^kp k ~ p c + VPsp'(ps)- 
k=l 

The super-critical sum is approximated by the integral, 

kr(k, t) dk. 
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Here, the lower limit is k = and not k c , because the characteristic scale, k, of k in r(k,t) 
is much larger than k c . Thus, the conservation identity takes the asymptotic form 



Notice that if r = 0, which corresponds to negligible super-critical clusters, (I6UJ) reduces to 
()54|) with 7] = 77*. Hence, r](t) has the effective initial condition 77(0) = 77*. 

In summary, the signaling problem for r(k,t) consists of the advection PDE (jHHj) . the 
creation BC 1)5 9 J) , and the functional dependence of the super-saturation, rj on r(k,t) in 
()60|). This signaling problem is nonlinear because rj, as a functional of r, appears in the 
advection PDE and, in addition, in the exponential creation rate in the BC. 

VI. CONCLUSIONS 

The discrepancy between the two accepted model of nucleation — the surface reaction 
model derived by Becker and Doring, and the diffusion limited growth model due to Lif- 
shitz and Slyozov — has now been resolved. Our new model, in which clusters interact with 
diffusing monomers by a surface reaction, predicts both models as limit cases. Although it 
provides a growth rate for clusters of all sizes, from small to large, its limiting behaviors are 
of special interest for us. In the limit of small clusters, the BD model emerges, and we can 
derive the Zeldovich creation rate of super-critical clusters using the BD kinetics. In the 
limit of large clusters, the diffusion limited growth model emerges, from which we derive the 
same evolution PDE for the density function, r(k,t), as LS. 

Not all our findings corroborate the classical assumptions and results. We find that the 
'common wisdom' about the monomer density at the cluster surface is wrong. In the classical 
DLG model, the surface density is chosen so that the surface reaction is in equilibrium. In 
our model, the balance between the surface reaction and diffusion leads to a surface monomer 
density substantially different from the 'equilibrium' value. Despite this, the new model still 
predicts the same growth rate as the classical DLG model (for large clusters). 

Another piece of 'common wisdom' is that BD pertains to small clusters and DLG to 
large ones. In our model, there is a new characteristic cluster size, k*, that separates small 
(k <C k*) from large {k 3> k*). This is useful in providing definite predictions for the validity 
of the model. For example, since nucleation happens around the critical size k c , and the 




(60) 
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Zeldovich creation rate assumes the small cluster limit of the model, the validity of the 
Zeldovich formula requires k c <C fc* . The case where k c is of the same order or larger than fc* 
is not covered by the current paper nor by classical Zeldovich nucleation theory. A separate 
investigation is required for this case, one that goes beyond or treatment of the monomer 
bath as a smooth 'mean field' and really accounts for its discrete and fluctuation nature. 

The monomer density around the cluster satisfies the diffusion equation and a mixed 
boundary condition at the surface of the cluster. To calculate the flux of the monomers 
into or out of the cluster, our model assumes that the monomer density in the vicinity of 
the cluster is a quasi-static solution of the diffusion equation. We show that the quasi- 
static 'surface' value of monomer concentration in our model is stable: starting from general 
initial conditions, the concentration of the monomer bath rapidly relaxes to the quasi-static 
approximation. 

The preceding insights are codified in a signaling problem for the cluster size distribution. 
It consists of an advection PDE consistent with diffusion limited growth, and a BC consistent 
with the Zeldovich creation rate. The PDE and BC contain the super-saturation as a 
parameter. The super-saturation, in turn, is a function of the cluster size distribution 
as dictated be the overall conservation of particles. Hence, the full signaling problem is 
non-linear. In particular, the Zeldovich creation rate is exponentially small as the super- 
saturation goes to zero. Therefore, accounting for small corrections to the monomer density 
is important, especially in the initial stage when new clusters are being created. 

Given an initial condition of pure monomers, quasi-static densities are formed during the 
ignition phase [3], whereby small, sub-critical clusters are created from the monomers. The 
ignition transient is a precursor to the nucleation of super-critical clusters. Therefore, once 
the nucleation 'has ignited', the monomers density is lower than the its original value, prior 
to the ignition transient. The signaling problem starts from an effective IC that explicitly 
accounts for the loss of monomers during ignition. This is especially relevant in our case, 
where the monomer-density is close to the saturation density, p s , and the dependence of the 
PDE and BC on the monomer density is very strong. 

Such is the state of the theory. We now conclude the conclusion by a small trespass 
into the domain of accountability: do the material parameters of real aggregation processes 
cooperate with the various assumptions of the modeling? We start with the characteristic 
cluster size fc*, most conspicuously present in our model. We examine physical constants 
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associated with an aqueous solution of calcium carbonate, CaC03. Its solubility is rather 
small, so the basic requirement of 'diluteness' is satisfied. Admittedly, CaCOs dissolves into 
positive (Ca 2+ ) and negative (CO 2- ) ions, which is not reflected in our aggregation model 
with identical particles. But we are examining crude order of magnitude estimates, so will 
not be distracted by the inconvenience, and we formally consider a positive-negative ion pair 
as a 'monomer.' 

Formula (f2T)j) for fc* contains molecular volume v, the saturation density of monomer p s , 
the diffusion coefficient D of monomer in solution, and the dissociation rate constant uj. 
E stlm a t e S of „ and D are read lly fonnd in a cnemrca! h -„„ t Q. We Us e „ a S an 
order of magnitude estimate of p s . The dissociation rate, u, is much more elusive. Here, 
we indulge in the activation energy model similar to Kelton's, in his review of glass-crystal 
transitions j^j]. The model is summarized by the formula 

D 

~2 



" = ^ e ""' ( 61 ) 



2/3 

Here, ^5-is the characteristic time for a monomer to diffuse a distance comparable to its 
own size, and (3 is an 'activation energy of dissociation', in units of ksT. Inserting (|61|) for 
us into (12*61). we find 



K = (3 • 16tt 2 ) (ePp s ) 



3 

) 



or, using p s = e , 

K = (3-167r 2 )e 3(/3 - a) . (62) 

Formula has other applications for us, besides the estimate (J62|) of fc*. For instance, 
recall that the quasi-static limit of the FBP (f2*Tj) is based on the 'diffusion time' much 
smaller than the dissociation time ^, so we required J^ ri 1. From ()61|) we get 

e' 3 (63) 



ojv 2 ' 3 

and so an activation energy ksT^ on the order of a few KbT is sufficient. 

Next, we examine the criterion |^ -C 1 for the validity of the Zeldovich nucleation rate. 
Using (fTUl I62J) . we find this implies a bound on the super-saturation, 

7] > cre Q_/3 . 
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Using the crude 'cube' model of bonding energy (JTJ) we estimate a by, cr « |a so the criterion 
on the supersaturation becomes 

T)^ae a - p . (64) 

A lower bound on the super-saturation might seem like a problem, as we expect rj to 
asymptote to zero in the long-time limit of an aggregation process. However, significant 
nucleation occurs at an early and relatively brief phase, so (JMj) should apply for the initial 
super-saturation. In late stage coarsening, 77 is much smaller than the RHS of (j64jl . but 
significant nucleation is not happening then. 

Numerical evaluation of fc*, or the RHS of (|64jl require actual values of a and ft. It is 
relatively easy to find a, as we have a = — log p s and p s is essenti ally the volume fraction 



of monomer in saturated solution. For instance, from the data in [14|, we find that 2.44 x 
10~ 4 cm 3 of solid CaC03 is soluble in 100cm 3 of water at room temperature. The volume of 
the solution is nearly 100cm 3 so the volume fraction occupied by Ca 2+ , CO 2 " is (roughly) 
p s = 2.44 x 10 -6 , and our estimate of a is a ~ 12.9. The bad news is that we do not know 
ft any better than we know uo in the first place, so we cannot do primia-facie evaluations of 
K nor the RHS of fU). 

Our policy is to use ()64j) to obtain bounds on ft for which our model is valid. For instance, 
we have seen that the quasi-static approximation requires ,9 > 1. In addition, our whole 
analysis is based on small super-saturation, so r/ <C 1. This is compatible with (fMj) only if 

ae a ~ f5 < 1. (65) 

Starting with a = 12.9, we find the LHS is unity if ft — a = 2.56, but if we increase the 
activation energy, ft, by 3ksT we get ft — a = 5.56 and the LHS of (|63|) is 0.061 <C 1. 
Inserting the value of ft — a = 5.56 into (JB2J) for k* results in k* w 8.3 x 10 9 . 

Our estimate of /c* does not really use values of D and v because the combination 
which appears in the formula ()26|) is absorbed by the Kelton formula (JoTj) for w. Again 
we 'turn the equations around': using physical values of D, v, and p s , and our 'made up' 
estimate of k*, we work backwards to w. From the Handbook of Chemistry and Physics 
[l^ . we find that the diffusion coefficients of Ca 2+ and CO 2 " are .792 x 10" 5 cm 2 /s and 
.9.23 x 10 _5 cm 2 /s. So, as a crude estimate we put D f=s 10 _5 cm 2 /s. The nucleous colume of 
solid CaCOs is v = 4.51 x 10~ 26 cm 3 . Inserting these values of D and v into (|26|). along with 
the previously determined p s and fc* we obtain the estimate of w, w ~ 2 x 10 6 s _1 . 
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APPENDIX A: APPENDIX 

1. Deterministic Cluster Growth According to BD 

Here we determine the domain of k in which the average change in the cluster size (8k) 
is a good approximation to the actual change in cluster size. BD assumes that the cluster 
grows and shrinks by means of discrete, independent Poisson processes. The cluster gains 
monomers at a rate c k p\ via the adsorption process, and via the emission process it loses 
monomers with rate d k . Therefore, the average change in size over a short interval St is 



Since the two controlling process are assumed to be independent Poisson processes, the 
mean square deviation of k from the average is 



The process can be considered deterministic if the deviation is much smaller than the ex- 
pected change, so 



(Sk) = (k(t + St) - k(t)) = (c kPl - d k )St, 



provided that St is small enough so that 



CfcPi — d k \ St <C k. 



(Al) 



{{Sk-(Sk)) 2 ) = {c kPl + d k )St 




(A2) 



Combining inequalities (|A1|) and ()A2|) we get a condition on St 



c^i + d k < \ (c k pi - 4)| St < \c k pi -d k \k. 



We can find St which satisfies this condition if 



CfcPi + d k < \c k pi - d k \k. 



(A3) 
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This does not hold for every value of k. In particular, we notice that for k = k c the RHS 
is zero. Hence, the k{t) that can be approximated as deterministic must be bounded away 
from k c . 

For k ^> 1 in (jA3|) we insert the approximations (0 EJ) for dk and and then examine 
its asymptotic form as rj — > + . The result is 



2 + 77 1 



<C rj 



k. 



(A4) 



For clusters with k ^> k c , this condition is satisfied trivially for rj <C 1. For k on the order 
of k c , we need to dig a little deeper. 

Since we are looking at the rj — > limit, we look at the leading order terms of (|A4)1 



2 

r]k 

For k near /c c we expand the RHS around k = k c . Keeping the leading order terms of the 
Taylor series of the RHS, we find an explicit condition for the evolution of clusters to be 
treated deterministically: k must be bounded away from k c according to 

- < \k - k c . 
V 

As rj — > 0, the excluded domain grows. On the other hand, its size relative to k c ~ 
shrinks to zero: 



677 



a-' 



Therefore, we see that for small super-saturation the evolution of clusters can be treated 
deterministically, for all but a vanishingly small domain around k c . 



2. Effective Super-Saturation, rj, After Ignition 

It has previously been shown that an initial condition of pure monomer leads to a transient 
'ignition' phase in which quasi-static densities of sub-critical clusters are created Q|. Since 
these sub-critical clusters are made of monomers, after the ignition transient the monomer 
density will be lower than the original value, prior to the ignition. We compute the value 
77* of super-saturation after the formation of the k < k c quasi-static densities, but before 
significant depletion by super-critical clusters. 
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The starting point if the sum ()39|) truncated to k < k c because the contribution from 
super-critical clusters is insignificant: 

p = ^kp k . 

l<k<k c 

Since the sub-critical cluster densities are quasi-static, we approximate them by the equilib- 
rium densities f!46|) . and the corresponding approximation of the sum is 

l<k<k c 

Recalling the definition of the super-saturation, (jHJ), pi = p s (l + ?7*), and taking the two-term 
expansion as 77* — > + , we find 

P~ Pc + V*J2k 2 p k s e~^ T . 

l<k<k c 

The sum on the RHS can be obtained by differentiating the series representation ()47j) of 
p(pi), setting pi = p s , multiplying by p s , and then truncating to k < k c . The derivative 
series converges at p\ — p s , so for 7] — > (and hence for k c — *> 00), the sum is asymptotic to 
PsP'(Ps)- Hence, 

P ~ Pc + V*Psp'(ps), 

and so, 

r] = 1 p ~ pc 
P'(ps) Ps 
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